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Abstract. - By combining exact diagonalization with scaling method, we show that it is 
possible to realize two channel spin | Kondo (2CK) conductance in a quantum dot at Coulomb 
Blockade, with an odd number of electrons and with contacts in a pillar configuration, as an 
applied orthogonal magnetic field B is tuned at an appropriate level crossing. 



A Quantum Dot (QD) weakly coupled to the contacts and tuned in a valley between two 
conduction peaks (Coulomb Blockade (CB)), is insulating if its charging energy is larger than 
the thermal energy However, when the number of electrons at the dot, N, is odd, below 
a characteristic temperature scale Tr-, a strongly correlated state between the dot and the 
contacts sets in, and the typical Kondo resonance in the conduction electron spectral density 
builds up at the chemical potential of the contacts fi. The striking result is that the linear 
conductance increases in the CB valley, when the temperature T is lowered, up to the unitarity 
limit 2e 2 /h for T = If N is even, the ground state (GS) of the QD is usually a spin 

singlet and the ordinary Kondo effect cannot occur. Nevertheless, it has been shown that level 
crossing between states at different S induced by a magnetic field B orthogonal to the dot 
can restore the degeneracy required for the Kondo effect to take place The occurrence 

of Kondo physics in quantum dots was predicted long ago in analogy to magnetic impurities 
in diluted metal alloys at very low temperatures R|. A dot at CB acts as a single magnetic 
impurity, but under controlled experimental conditions. 

A realistic description of the Kondo effect in alloys has to take into account that the 
impurity states carry also orbital angular momentum together with spin momentum. Hence, 
electrons with different orbital momentum can access the impurity and both multi-orbital 
1-channel Kondo Effect (1CK) and many-channel Kondo effect (MCK) can take place. Total 
angular momentum is conserved in the scattering, as it is appropriate for an atomic impurity 
in an isotropic metal ||. 

Similarly to the ordinary 1CK, MCK shows a logarithmic low-temperature raise in the 
resistivity before the perturbative expansion breaks down (T ~ Tr-). However, if the number 
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Fig. 1 - (a): Sketch of the dot. (b) Sketch of the conductance as a function of the temperature T and 
the gate voltage V g . At 2CK non Fermi liquid point, the conductance shows a y/T dependence. 



of channels k > 2S (overscreening), the physics of the two systems at T <C Tk is dramatically 
different. The overscreened impurity interacts again antiferromagnetically with the other 
conduction electrons. Accordingly, the strongly coupled antiferromagnetic (AF) fixed point 
becomes unstable at T — and moves to intermediate coupling. The Fermi-liquid phase 
breaks down and leading corrections to the resistivity p(T) are characterized by a non-Fermi 
liquid square-root temperature dependence H . The prototype of this situation is an S = | 
impurity interacting with two channels of conduction electrons. 

Possibly, 2CK zero-bias conductance anomaly has been detected in transport measure- 
ments in Ti— and Cu— nanoconstrictions with controlled doping, with Tk between 0.5 K 
and 4.6 K jL0|. Also, interstitial atoms in glassy metals could be described as tunneling 
two-level systems (TLS) and they could give rise to 2CK by scattering with conduction elec- 
trons 11 1 . However, a careful estimate of the Kondo temperature in these materials is still a 
widely debated question and appears to give a value below 10~ 4 K 

In this letter we propose to search for 2CK spin | in a quantum dot (QD) biased at 
Coulomb blockade with odd N and total spin S = 1/2. We discuss the geometry, the working 
conditions, and the most appropriate energy level structure to achieve 2CK within the limits 
of present technology. 

We consider a QD in a pillar configuration with cylindrical symmetry about the z-axis (see 
Fig.l a) ). We label dot states with total energy N Ei ( i = for the GS) and total angular 
momentum M along the z axis. By means of single-electron tunneling with the contacts, 
the dot can exchange orbital momentum m and/or spin a with conduction electrons. We 
describe a particular scenario derived from results of an exact diagonalization of electrons in a 
two-dimensional parabolic confining potential with Coulomb interaction energy scale U |l3| ]. 
The addition energies of the dot can be shifted with respect to p by tuning a gate voltage V g . 
The field B orthogonal to the dot can induce strong orbital effects, which may eventually give 
raise to crossings between states with different M and S. The dot is tuned at CB with N = 5 
at B* where two S = \ levels with M = 4 and M = 6 cross (Fig.2 ). 
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Fig. 2 - The lowest dot energy levels Ei, obtained by exact diagonalization, are plotted as a function 
of the magnetic field B. At the degeneracy point B* ~ 5 Tesla two S = ^ levels for N=5 cross. The 
dot switches between the two states via electron cotunneling to and from the contacts in the particle 
(p) and the hole (h) channel ( marked by thick arrows: the starting and ending points are displaced 
from B* for clarity). Intermediate virtual many-body states with N=4,6 are singlets. Symbols are 
explained in the text. 

Within the subspace S spanned by the states \N, M,S Z ) = |5,4, ±|) and |5,6,±|), the 
QD may scatter a conduction electron with spin er from m = to m = 2, and vice versa. In 
this respect, the dot acts as an impurity spin S = 3/2 in interaction with two channels of 
conduction electrons (In our case the channel index is the spin index a =f,J,). This would 
attempt to realize underscreening of the dot's spin with Kondo temperature T C ross- However, 
if Zeeman spin splitting (ZSS) Az > T CIOSS (fc^ = h = 1 throughout the paper), two of the 
dot states are ruled out at temperatures T such that Az > T > T cross . The system may flow 
toward an overscreened S = | 2CK fixed point, provided its Kondo temperature T 2 ck is such 
that T 2 ck > 2c ross . (The ZSS should not affect the density of conduction electron states of 
both spin polarization at fi.) 

We now discuss for which values of the parameters the system may exhibit 2CK and 
estimate the corresponding Kondo temperature (see Fig. 3) by combining exact diagonal- 
ization results for an isolated dot with N = 4,5,6 electrons with scaling approach. We 
construct the even-parity single-particle states labeled by qma out of the wavefunctions 
e «(fcF,i m +g)2<^ e ^ m ^(p) [Q ; delocalized in the z-direction, with e ~ \x and q <C k F ^ m (k F ^ m is 
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the Fermi wavevector of the 1-d subbands of motion along z, labeled by I and linearized about 
//). We denote the corresponding operators by b qma , b qma and drop the energy and subband 
index for simplicity. The Hamiltonian for the leads is given by H b = '}2 qma {^qm~ ft)b qma b qma . 

Dot states can be expanded in the basis of Slater determinants of single particle orbitals 
(n, m ) associated to the fermion operators d nm<7 , dj imcr (in a in-plane parabolic confining po- 
tential, n = 0, 1, ..oo and m = —n, —n + 2, ...,n—2,n). In Fig. 2 we represent the determinants 
with S z = \ giving the largest contribution to the N — 4,5,6 GS, by filling a sequence of 
boxes, n increases by one along the vertical rows of boxes , while m increases along the hori- 
zontal rows. A large magnetic field B || z favors the occupancy of positive m states. All the 
states have the lowest possible S, since Hund's rule does not apply in the presence of a rather 
strong B. The larger is U, the higher in energy are N — 4, 6-triplet states. Therefore, we 
neglect virtual transitions to triplets, since B should be cranked up a lot in order to induce 
transitions to higher spin states of the dot. 

Within the subspace S, the Hamiltonian of the isolated dot is: 

H d = Kj2s z aa + A z (P n -P u ) (1) 

a 

In Eq.(0), we have introduced the following dot operators: 

P aa ' = 5 £ MOW I ; S% = ±5>a<7)r° Q ,<5a' a \ 

a aa 

a labels states with different M and r^ a , denote the Pauli matrices with a = 1, 2, 3. The first 
term of eq.([l]) is a possible detuning of B from B* , which breaks the degeneracy of the dot 
states at different M and is ignored in the following. The second term is the ZSS term. 

The model hamiltonian for the tunneling between dot and leads is Ht — T J2 nqm a ^nma-^qmc7+ 
h.c. To lowest order in |r| 2 , H t connects dot states S 3 to states with N ± 1. At CB, 
6 Eo — 5 Eq < ii < 5 Eq — 4 Eo and the two most relevant cotunneling processes giving raise to 
tunneling of the dot between states with M = 4 and M = 6 are marked by thick solid arrows 
in Fig. 2. The labels p, h refer to whether the lead intermediate state has one extra electron 
(particle), or one extra hole (hole). 

Let us denote the states of the contacts by \[Af], 4>), where Af is the number of electrons and 
<f> is a collective label for all the appropriate quantum numbers. Generalized magnetization 
vector operators acting on the contacts at the dot site z = in the perturbative h, p-process Jl5| 
are given by: 



er (7 



qmm 



[Af±l T {qma}],<f>)\[M±l],<j>' 
[M±{qm'a}},A([W,4," 



(2) 



Generalized scalar charge density operators, fP' , , are obtained by substituting T mm ' with S mm > 
in Eq.(||). A Schrieffer- Wolff transformation provides the effective Hamiltonian restricted to 
the subspace S: 
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+ § aa ■ m h ss \ + y J2J2 § °*-™ X ™ ; ( 3 ) 

x—p.h cr 

where a = —a and J ~ |r|7| 5 £ - 4 £ - J ~ |r| 2 /| 5 ^ - 6 B + A*]. 

The dominant p and h processes are diagonal in the spin indices and involve only lead 
electrons with a opposite or equal to S z , respectively (see Fig. 2). The other possibility 
corresponds to subleading processes in which the intermediate dot state is a triplet. Y is the 
coupling strength for tunneling processes that change both M and S z (off-diagonal terms in 
the spin labels). 

We calculate the couplings in H^s by means of exact numerical diagonalization (Inset of 
Fig.(||) ). The transverse and longitudinal amplitudes, ( J z , J ± ), ( J z , J ± ), (Y z , Y^), are all 
of the same order of magnitude. Anisotropies between transverse and longitudinal couplings 
are irrelevant to our argument anyway. Therefore, we assume isotropic magnetic couplings in 

T > T cross > A z case : 

A scaling analysis on i?Eff is performed by rescaling the band cutoff as D — > D — 5D and 
by working out the corresponding renormalization of the coupling strengths J, J and Y. This 
provides third order scaling equations given by: 

{-v{0)[J 2 + Y 2 } + v 2 {0)[J 3 + JY 2 }} ■ 
{-v{0)[J 2 + Y 2 } +v 2 {0)[J 3 + JY 2 }} ■ 

-- {-v(0)[JY + JY] + v 2 {0)[J 2 + J 2 ]Y} (4) 

(v(0) is the density of states at \x for each spin polarization). 

The initial values of the parameters Jo, Jo, Yq define a crossover temperature to strong 
coupled regime, T cross . Because of the spin flip term proportional to Y in Eq.(||), diagonal- 
ization of Hes selects one preferred scattering channel, which does not conserve spin and 
orbital momentum separately. Hence, below T cross , the system flows towards a 1-channel 
underscreened strongly-coupled state. Given isotropic couplings, Jo ~ Jo ~ Yo, T CIOSS = 

D oy /v{p)J exp[--s$p 5 ]- 
Az > r cr o SS case : 

Instead, we are interested in the complementary situation, Az > T > T cross . The two 
S z — — \ dot's states decouple from the scaling when T < Az, and spin-i 2CK may take place 
in the system. At temperatures below the decoupling point H c g no longer applies. Indeed, 
the off-diagonal term proportional to Y in Eq.(^) does not contribute to scaling anymore. As 
the dot lays within either state with S z = h, becomes a spin ^-operator, — i, while 
Pll = 0. The effective Hamiltonian describing low-energy dynamics at the scale Az takes 
the form of a 2CK-Hamiltonian Hick-, with the spin index playing the role of the channel 
index. Coupling constants in Hick are given by the coupling strengths J, J scaled down to 
T = Az, J*, J*. After a particle- hole transformation is performed in the j-spin channel only 
{bqm] ( — l) m ^mt' w ^ n TR = Q,2 and m m), the Kondo interaction term in H2CK is 
unchanged, while a sign is reversed in the potential part, yielding: 



d.J 



dJ 



d Hw ) ' 

dY 



H 2 ck = 2( J VtT + J * Pll) + J *S ■ "T-TT + J *S ■ mil 



(5) 
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As shown in the inset of Fig. 3, fine-tuning of dot's parameters is possible, to achieve the 
symmetry condition J* = J*. Due to the transformation all the virtual processes are all of the 
same kind and the label p, h has been dropped. The potential scattering in Eq.(^J) generates 
a phase shift responsible for depletion of p states at the Fermi surface. For this reason, it 
has to be taken into account by diagonalizing it at each step of the renormalization process 
described by the flow equation: 

d J* 

-^V-=-^(0)(J*) 2 + K(0)) 2 (J*) 3 , (6) 

(f p (0) is the p density at the Fermi level) jL6). From the scaling equation Eq.(||), one derives 
the corresponding 2CK-temperature: T 2 qk = Az^ p (0)J* e~ " p (o>j* . Below, we show that 
there is a wide range of parameters for which the condition T 2 qk > T cmss is realized. In this 
case, as T — > 0, the system may flow toward a "spin-1/2" 2CK fixed point. Indeed, spin | 
2CK takes place if coupling of the two conduction channels is symmetric, J* — J* . Although 
this condition is quite demanding, it should be possible to achieve it by properly tuning the 
gate voltage V g , as shown in the inset of Fig. 3. Inclusion of the excited states we neglected 
in our work (N = 4,6 triplets and so on), will shift the value of V g corresponding to the 
symmetry point. The signature of 2CK spin i fixed point will be in a \T dependence of the 
conductance below T 2 qk at a special V g point, instead of the T 2 behavior, typical of a Fermi 
liquid fixed point |^](see Fig.l b)). In Fig. 3 we plot the corresponding relevant temperature 
scales, T 2 ck and T cross vs. the hybridization parameter between dot and leads, A = i/ p (0)|r| 2 , 
in units of the single particle level spacing lu = \/lu^ + uj 2 /A =6 mcV (uid = AraeV is the 
frequency associated to the confining potential and lo c = 9meV is the cyclotron frequency at 
B = B*). Appropriate choice of the barrier thickness can fix the hybridization parameter A 
within the required range Az > T 2 qk > T CIOSS . 

In conclusion, we propose to search for orbital "spin i" 2CK in a QD in a pillar structure, 
tuned at CB with an odd number of electrons and an appropriate magnetic field along the 
axis, corresponding to the degeneracy point between GS levels with different values of M. 
Detecting the effect requires an appropriate control of the hybridization of the dot with the 
contacts and a proper tuning of the gate voltage V g . A measurement of the conductance G(T) 
vs. T at low temperature will exhibit quadratic behavior and will cross over to a square root 
behavior, as V g is tuned so that the two channels become symmetric. 

* * * 
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tion parameter A, according to expressions given in the text. The 2CK conductance may occur when 
Az > Tick > T CIOSS (hatched area). Values of the parameters are: uio =6 meV, £>o=30meV, Az — 
1.4 K, v(0) =0.084meV _1 . Inset: Plot of the couplings J z and J z when the position of the dot levels 
relative to n is shifted by changing V g . 
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